Risk Analysis, Vol. 30, No. 7, 2010 


DOI: 10.1111/j. 1539-6924.2010.01427.x 


Development of a Dose-Response Model 
for SARS Corona virus 
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In order to develop a dose-response model for SARS coronavirus (SARS-CoV), the pooled 
data sets for infection of transgenic mice susceptible to SARS-CoV and infection of mice 
with murine hepatitis virus strain 1, which may be a clinically relevant model of SARS, were 
fit to beta-Poisson and exponential models with the maximum likelihood method. The ex¬ 
ponential model (k = 4.1 x 10 2 ) could describe the dose-response relationship of the pooled 
data sets. The beta-Poisson model did not provide a statistically significant improvement in 
fit. With the exponential model, the infectivity of SARS-CoV was calculated and compared 
with those of other coronaviruses. The does of SARS-CoV corresponding to 10% and 50% 
responses (illness) were estimated at 43 and 280 PFU, respectively. Its estimated infectivity 
was comparable to that of HCoV-229E, known as an agent of human common cold, and also 
similar to those of some animal coronaviruses belonging to the same genetic group. More¬ 
over, the exponential model was applied to the analysis of the epidemiological data of SARS 
outbreak that occurred at an apartment complex in Hong Kong in 2003. The estimated dose 
of SARS-CoV for apartment residents during the outbreak, which was back-calculated from 
the reported number of cases, ranged from 16 to 160 PFU/person, depending on the floor. 
The exponential model developed here is the sole dose-response model for SARS-CoV at 
the present and would enable us to understand the possibility for reemergence of SARS. 


KEY WORDS: Dose-response model; maximum likelihood method; SARS coronavirus; SARS out¬ 
break in Hong Kong 


1. INTRODUCTION 

The causal pathogen of severe acute respira¬ 
tory syndrome (SARS) is a newly isolated coron¬ 
avirus (SARS-CoV) that first appeared in late 2002 
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in Guangdong Province, People’s Republic of China. 
In the spring of 2003, a large outbreak of this se¬ 
vere pneumonia occurred in Hong Kong and rapidly 
spread throughout the world. Ultimately, 8,096 cases 
of SARS were identified in 29 countries or areas and 
774 patients reportedly died.* 1 ’ The rapid transmis¬ 
sion and high mortality rate made SARS a global 
threat for which no efficacious therapy was available 
and empirical strategies had to be used to treat the 
patients. 

Since the SARS pandemic in 2003, develop¬ 
ments of vaccines' 2 ’ and antivirals^ have been 
rapidly proceeding for treatment of SARS patients 
and prevention of its reemergence. In parallel, 
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mathematical models expressing the SARS propaga¬ 
tion from person to person have been developed to 
discuss countermeasures against such a disease trans¬ 
mission from community to national scales. (4,5) These 
models successfully explained the epidemic curve by 
estimating model parameters of infectious rate, re¬ 
moval rate, and so on. However, these models did not 
incorporate primary exposures or transmission via 
environmental reservoirs and the model parameters 
would not be consistent and highly influenced by nu¬ 
merous factors in environment and human behavior. 
Such models with parameters estimated on a case-by- 
case basis are of little use in any outbreaks occurring 
in other situations. On the other hand, dose-response 
models to characterize the interaction between hu¬ 
man and virus seems robust and applicable to as¬ 
sessing the risk of SARS via any possible routes of 
infection (e.g., aerosols containing virus particles, 
surface or hand contamination). Nevertheless, no 
dose-response model for SARS-CoV is available 
at the moment mainly due to unavailability of 
data sets challenging humans or animals with this 
virus. 

In this article, we develop a dose-response model 
for SARS-CoV based on two data sets for infection 
of transgenic mice susceptible to SARS-CoV and in¬ 
fection of mice with murine hepatitis virus strain 1 
(MHV-1) that may be a clinically relevant model of 
SARS/ 6 ’ And also, as an example of model applica¬ 
tion, we analyze the epidemiological data of a SARS 
outbreak that occurred at an apartment complex in 
Hong Kong in 2003 with the developed model. 

2. METHODS 

2.1. Data Sets for Model Development 

Coronaviruses cause acute and chronic respi¬ 
ratory, enteric, and central nervous system (CNS) 
diseases in many species of animals including hu¬ 
mans. Previous to the emergence of SARS, there 
were two prototype human coronavirues OC43 and 
229E, both etiologic agents of the common cold. 
SARS was the first example of serious illness in 
humans caused by coronaviruses/ 7 ’ Since the out¬ 
break in 2003, many researchers have worked to 
elucidate structures of viral genes and proteins, 
mechanisms of infection and replication, and its 
pathogenesis. Coronaviruses are divided into three 
groups (groups I to III) based on the genome 
sequences and SARS-CoV may be a member of 
group II as well as murine hepatitis virus (MHV), 
bovine coronavirus, porcine hemagglutinating en¬ 


cephalomyelitis virus (HEV), equine coronavirus, 
and human coronavirues OC43 and NL63, which also 
cause respiratory infections. Among them, MHV 
that infects both mice and rats often has been studied 
as a suitable model of human coronavirus diseases. 

Coronaviruses are generally restricted in their 
host range and viruses associated with disease in one 
species can be limited in their ability to replicate in 
other species/ 8 ’ SARS-CoV differs from this gen¬ 
eral pattern. This virus infects and replicates in mice, 
ferrets, hamsters, cats, and several species of non¬ 
human primates (cynomolgus and rhesus macaques, 
African green monkeys, and marmosets). Neverthe¬ 
less, most attempts to reproduce completely human 
clinical disease and pathological findings in these an¬ 
imals failed. On the other hand, De Albuquerque 
et al. l9> demonstrated that intranasal infection of 
A/J mice with MHV-1 produced pulmonary patho¬ 
logical features of SARS. From the fact that all 
MHV-l-infected A/J mice developed progressive in¬ 
terstitial pneumonitis, including dense macrophage 
infiltrates, giant cells, and hyaline membranes, result¬ 
ing in death of all animals, they concluded that A/J 
mice infected with MHV-1 would be a potentially 
useful small animal model of human SARS that de¬ 
fines its pathogenesis and suggests treatment strate¬ 
gies. We employed the data set (2 in Table I) chal¬ 
lenging A/J mice with MHV-1 as a surrogate data set 
for SARS-CoV. The data set was obtained by mon¬ 
itoring the survival of four groups of the mice (n = 
5 per group) for 21 days after intranasal inoculation 
with MHV-1 at 5, 50, 500, and 5,000 PFU, respec¬ 
tively. In other approaches, transgenic (tg) mice ex¬ 
pressing the human receptor for SARS-CoV, which 
are very susceptible to SARS-CoV, have been devel¬ 
oped and used for pathogenesis studies/ 10 ’ DeDiego 
et a// 11 ’ reported a data set for infection of the tg 
mice with recombinant SARS-CoV (1 in Table I). 
In their experiment, four groups of the tg mice (n = 
2 to 6 per group) were intranasally inoculated with 
240, 800, 2,400, and 12,000 PFU of rSARS-CoV, re¬ 
spectively, and the survival was monitored 13 days. 
The details of data sets 1 and 2 are shown in the 
Appendix. 

The endpoint was death of mice for both data 
sets 1 and 2. Since neither infection nor illness of mice 
has been analyzed in these experiments, we assumed 
that all mice with the illness died. In general, mortal¬ 
ity of animals with severe illness like SARS would 
be dependent on their physical strength, which in¬ 
dicates a potential of recovery from the illness, and 
availability of medical treatment. As described 
above, a high mortality was observed during the 
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Table I. Data Sets on Dose-Response Relationship for Coronavirus Infection via Intranasal Route* 


No. 

Virus 

Genetic 

Group 

Host 

Reported 

Symptom 

Number 

of Doses 

Endpoint of 
Response 

Reference 

Animal 

Age 

1 

rSARS-CoV** 

2 

tgMicetf 

? 

R.CNS 

4(1ft) 

Death 

11 

2 

MHV-1 

2 

Mice 

6 to 8 weeks 

R 

4(2) 

Death 

9 

3 

HCoV-229E 

1 

Humans 

18 to 50 years 

R 

4(4) 

Illness (Cold) 

12 

4 

MHV-S 

2 

Mice 

3 days 

CNS.H 

4(1) 

Death 

13 

5 

MHV-S 

2 

Mice 

1 week 

CNS.H 

4(1) 

Death 

13 

6 

MHV-S 

2 

Mice 

2 weeks 

CNS.H 

4(1) 

Death 

13 

7 

MHV-2 

2 

Mice 

3 weeks 

H 

3(0) 

Death 

14 

8 

MHV-2 

2 

Mice 

4 weeks 

H 

3(0) 

Death 

14 

9 

HEV-67N 

2 

Mice 

1 week 

CNS 

4(1) 

Death 

15 

10 

HEV-67N 

2 

Mice 

4 weeks 

CNS 

5(1) 

Death 

15 

11 

HEV-67N 

2 

Mice 

8 weeks 

CNS 

5(1) 

Death 

15 

12 

HEV-67N 

2 

Rats 

1 week 

CNS 

4(1) 

Death 

16 

13 

HEV-67N 

2 

Rats 

4 weeks 

CNS 

3(1) 

Death 

16 

14 

HEV-67N 

2 

Rats 

8 weeks 

CNS 

3(2) 

Death 

16 

15 

IBVA-5968 

3 

Chicks 

9 weeks 

R,K 

6(4) 

Death 

17 


R = respiratory; CNS = central nervous system; H = hepatitis; K = kidney damage. 
*Only data set 15 was obtained by intratracheal inoculation. 

“Recombinant SARS-CoV. 

^Transgenic mice expressing the SARS-CoV receptor. 

^Number of dose points corresponding to other than 0 or 100% response. 


SARS outbreak, although the patients were treated 
with some empirical strategies. In the experiments 
using mice, higher mortality can be expected since 
the mice developing illness had never received 
any medical treatments and basically the physical 
strength of mice would be much lower than that of 
humans. This could support the above assumption 
and therefore we decided to use the data sets 1 and 
2 to develop the dose-response model with the end¬ 
point of human illness, that is, SARS. 

In order to compare the infectivity of SARS- 
CoV with other coronaviruses, we also collected data 
sets (3 to 15 in Table I) challenging humans and an¬ 
imals with HCoV-229E, MHV, HEV, and infectious 
bronchitis virus (IBV), which had been studied for 
many years before the emergence of SARS. Except 
for one data set (15 via tracheal inoculation), all these 
data sets were obtained by the intranasal inoculation, 
which is probably the primary route of infection with 
SARS-CoV. 


2.2. Fitting of Dose-Response Models to Data Sets 

Each data set in Table I was fitted to two types 
of dose-response model, that is, exponential and 
beta-Poisson models/ 18 ) with maximum likelihood 
method and goodness-of-fit to the models compared 


based on their likelihoods. The general forms of 
these models are as follows: 


Exponential model: p(d) = 1 — exp ( — — 


Beta-Poisson model: p(d ) = 1 — 


1+ aH 24 *- 1 

A50 


where p{d) is the risk of illness at the dose of d\ and C, 
V 50 , and a are parameters specific for the pathogen. 
Parameter k in the exponential model equals the 
reciprocal of the probability that a single pathogen 
will initiate the response. Parameters V 50 and a in 
the beta-Poisson model are the median dose to get 
the response and the exponential fitting parameter, 
respectively. 

The beta-Poisson model was fitted to only four 
data sets 2, 3, 14, and 15 since it is not meaningful 
to fit this model to data sets including less than two 
dose points corresponding to other than 0 or 100 % 
response. The model fitting was done via maximum 
likelihood estimation, using a quasi-Newton method, 
implemented in the package of R version 2.6.2 for 
Windows/ 19 * The sensitivity analysis of parameters 
of the best-fitted model was performed with 10,000 
bootstrap trials by a program written on R to esti¬ 
mate the doses corresponding to 10% and 50% re¬ 
sponses (ID 10 and ID 50 ). 






1132 


Watanabe et al. 


2.3. Model Application: SARS Outbreak at an 

Apartment Building in Hong Kong 

During the pandemic of SARS in 2003, Hong 
Kong was the hardest hit reporting area with 1,755 
cases and 299 deaths in a population of 6.7 million.* 
The large community outbreak at an apartment com¬ 
plex named Amoy Gardens affected more than 300 
among 20,000 residents in the early stage of the out¬ 
break in Hong Kong. According to the report by 
Department of Health (DOH), Hong Kong SAR, 
the outbreak was begun with an index case who vis¬ 
ited his family living on the 16th floor in Block E of 
Amoy Gardens.* 20 ' 21 ) He stayed overnight and used 
the toilet there. DOH concluded that the sewage con¬ 
tamination associated with the index case (and other 
infected persons) excreting coronavirus that gained 
entry to households through the bathroom floor 
drain with dried U-traps was the primary cause of the 
SARS outbreak. The building’s ventilation system 
made a significant contribution to the entry of virus 
into households. According to McKinney et al. ,* 22 * 
many residents had installed high-powered fans in 
their small bathrooms with capacities 6 to 10 times 
higher than the required capacity and the fans cre¬ 
ated large negative pressure and drew air from waste 
pipes. The WHO environmental team verified that 
sewer gas and aerosolized droplets, which the hy¬ 
draulic action caused by flushing toilets generated, 
were being drawn into the bathrooms from the waste 
pipe system. Although DOH mentioned other pos¬ 
sible causes of the SARS outbreak, such as person- 
to-person transmission, vectors acting as mechanical 
carriers for the virus, and environmental contamina¬ 
tion, we assumed that all cases were infected via the 
airborne transmission, reported as the primary cause, 
in our model application. 

DOH investigated epidemiological data of 321 
confirmed SARS cases. Although the detailed data 
are not reported in the literature, the location of flats 
for 99 cases in Block E has been published in the 
modeling study by Li et alS 21 ' 1 Since the number of 
residents at each flat during the outbreak has not 
been investigated, their group assumed that each flat 
housed four persons, which was the largest number 
of cases in any one flat, and calculated the risk of ill¬ 
ness, equivalent to the attack rate, as the number of 
cases divided by the assumed number of residents.* 23 * 
By using the dose-response model developed here, 
we estimated the dose ( d ) of SARS-CoV for cases in 
Block E from the attack rate (p(d )). The dose esti¬ 
mated on the above assumption might be lower than 


the true value because number of residents in some 
flats would be below four and there were possibly un¬ 
occupied flats. We also calculated the attack rate as 
the number of flats having at least one case divided 
by the total number of flats and estimated the dose, 
which would be higher than the real value, from the 
attack rate. 

3. RESULTS 

3.1. Dose-Response Model for SARS-CoV 

Table II shows the result of fitting models to the 
data sets. The fits of the beta-Poisson model to data 
sets 2, 3, 14, and 15 could not be rejected on statis¬ 
tical grounds, as indicated by pm values greater than 
0.05. Similarly, the exponential model provided sta¬ 
tistically significant fits to all data sets except 14 and 
15. As for data sets 2 and 3, both beta-Poisson and 
exponential models were acceptable. We tested the 
null hypothesis that the beta-Poisson model provided 
a statistically significant improvement in fit to these 
data sets rather than a more parsimonious (expo¬ 
nential) model by comparing the difference of vari¬ 
ance between two models against the critical value 
(3.84) on the chi-square distribution with 1 degree of 
freedom. As the result, since the hypothesis was re¬ 
jected for both data sets, the exponential model was 
used for the analysis of viral infectivity in the latter 
part. 

The estimated parameter k in the exponential 
model for rSARS-CoV (data set 1) was very close 
to that for MHV-1 (data set 2). In order to decide 
whether these data sets could be pooled or not, the 
following statistic A was calculated and compared 
with the critical value (3.84) in the chi-square distri¬ 
bution with 1 degree of freedom:* 18 ) 

A = Y p -(Y 1 + Y 2 ), 

where Y i, Y 2 , and Y p are deviances of fits of expo¬ 
nential model to data sets 1 and 2 and pooled data 
sets, respectively. As shown in Table II, we yielded In 
k — 6.01 and Y p — 1.75 when the exponential model 
was fitted to the pooled data sets (data set p). The 
value of A (0.17) smaller than 3.84 indicated that data 
sets 1 and 2 could be pooled to p. Consequently, the 
exponential model (Fig. 1) with k — 4.1 x 10 2 was 
employed as the dose-response model for SARS due 
to intranasal infection. 
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Table II. Parameters and Likelihood of Dose-Response Models Fitted to Data Sets on Coronavirus Infection 


No. 

Virus 


Host 


Beta-Poisson(BP) 


Exponential (EXP) 

Better Fit 
Model** 

Unit of 

Dose 

Species 

Age 

In a 

In V 50 

Deviance 

Ait 

In k 

Deviance 

Ait 

1 

rSARS-CoV 

tgMice 

? 

_* 

- 

_ 

_ 

5.82 

0.97 

0.81 


PFU 

2 

MHV-1 

Mice 

6 to 8 weeks 

1.23 

5.71 

0.54 

0.76 

6.15 

0.61 

0.89 

EXP 

PFU 

Pt 

SARS-CoV 

Mice 

? 

6.50 

5.64 

1.75 

0.94 

6.01 

1.75 

0.97 

EXP 

PFU 

3 

HCoV-229E 

Humans 

18 to 50 years 

- 0.21 

2.38 

1.42 

0.49 

2.92 

2.31 

0.51 

EXP 

TCD 50 

4 

MHV-S 

Mice 

3 days 

- 

- 

- 

- 

6.32 

0.56 

0.91 


PFU 

5 

MHV-S 

Mice 

1 week 

- 

- 

- 

- 

8.62 

0.56 

0.91 


PFU 

6 

MHV-S 

Mice 

2 weeks 

- 

- 

- 

- 

13.6 

0.24 

0.97 


PFU 

7 

MHV-2 

Mice 

3 weeks 

(No death among 5 mice at doses from 7 

x 10 2 to 7 x 

10 4 ) 


PFU 

8 

MHV-2 

Mice 

4 weeks 

(No death among 5 mice at doses from 6 

x 10 3 to 6 x 

10 s ) 


PFU 

9 

HEV-67N 

Mice 

1 week 

- 

- 

- 

- 

2.39 

0 ++ 



PFU 

10 

HEV-67N 

Mice 

4 weeks 

- 

- 

- 

- 

7.62 

0.62 

0.96 


PFU 

11 

HEV-67N 

Mice 

8 weeks 

- 

- 

- 

- 

7.62 

0.62 

0.96 


PFU 

12 

HEV-67N 

Rats 

1 week 

- 

- 

- 

- 

3.08 

o+t 



PFU 

13 

HEV-67N 

Rats 

4 weeks 

- 

- 

- 

- 

8.66 

0.54 

0.76 


PFU 

14 

HEV-67N 

Rats 

8 weeks 

-1.04 

7.42 

1.97 

0.16 

9.34 

6.36 

0.04 

BP 

PFU 

15 

IBVA-5968 

Chicks 

9 weeks 

-5.11 

109 

2.93 

0.57 

11.4 

26.6 

< 0.01 

BP 

CD 50 


*Beta-Poisson model was not fitted to data sets including less than two dose points corresponding to other than 0 or 100% response. 
**Beta-Poisson model was employed if it provided a significant improvement of fit (p > 0.05) rather than exponential model. 

TData set p pooled data sets 1 and 2. 

^Zero deviance means that all data were on the model. 



Fig. 1 . Exponential dose-response model fitted to the pooled data 
sets for SARS-CoV. 

3.2. Infectivity of SARS-CoV 

Fig. 2 illustrates the doses (IDio and ID 50 ) of 
SARS-CoV and other coronaviruses corresponding 
to 10% and 50% responses calculated with the dose- 
response models listed in Table II. IDi 0 and ID 50 
of SARS-CoV were 43 (95% Cl = 20 to 81 PFU) 


and 280 PFU (95% Cl = 130 to 530 PFU), respec¬ 
tively. The doses (IDi 0 = 2.0 TCD 50 ; ID 50 = 13 
TCD 50 ) for HCoV-229E seem relatively low; how¬ 
ever, Schmidt et n/. (24 ) reported that the 50% end¬ 
point assay was about 10 to 30 times less sensi¬ 
tive than the plaque assay with HCoV-229E. Ac¬ 
cordingly, the estimated infectivity of SARS-CoV 
would be comparable to that of HCoV that causes 
a mild cold in humans. Similar infectivity was also 
observed in data sets of three-day-old mice infected 
with MHV-S (data set 4) and eight-week-old rats in¬ 
fected with HEV-67N (data set 14). These animal 
coronaviruses belong to the same genetic group as 
SARS-CoV. The infectivity of coronavirus might be 
related to the viral evolution, although a lot of un¬ 
known factors exist. The effect of host age was ob¬ 
served in data sets for MHV-S (data sets 4 to 6 ) 
and HEV-67N (data sets 9 to 14) and it was obvi¬ 
ous that young mice and rats are more susceptible 
than old ones. In contrast, infection among children 
and adolescents was relatively uncommon^ and the 
mortality of senior people was comparatively high in 
the SARS pandemic .^ This age-dependency of hu¬ 
man susceptibility to SARS-CoV is different not only 
from other coronaviruses but also from other human 
respiratory viruses. Avian influenza A (H5N1) virus 
is a highly virulent respiratory virus like SARS-CoV. 
The epidemiologic investigations show that this virus 























1134 


Watanabe et al. 



Fig. 2. Doses (IDio and ID 50 ) of SARS-CoV and other coron- 
aviruses corresponding to 10 and 50% response. Error bars indi¬ 
cate 95% confident intervals. 

primarily infects young people and Goicocchca* 27 ' 
made a comment that the receptor recognized by 
this virus, which is expressed in the lower respiratory 
tract, may be expressed in the upper airway in chil¬ 
dren, increasing the risk of infection. In case of res¬ 
piratory syncytial virus (SRV), another well-known 
human respiratory virus, the serious forms of the dis¬ 
ease (principally bronchiolitis and interstitial pneu¬ 
monia) are found most frequently in infants under six 
months of age and the disease is progressively milder 
with increasing age.* 28 ' 

3.3. Back-Calculation of Viral Dose for Residents 

at Amoy Gardens 

The doses of SARS-CoV for residents in Block 
E, Amoy Gardens, estimated based on the assumed 
number of residents and the number of affected flats, 
are described in Table III. As Yu et a/.* 23 ' analyzed, 
the attack rate was dependent on the floor and the 


middle floors between levels 14 and 23 had the high¬ 
est risk since the flat where the index case visited 
was located there. Therefore, the estimated dose of 
SARS-CoV (63 to 160 PFU) for residents on the mid¬ 
dle floor was also higher than the others. The es¬ 
timated dose (42 to 117 PFU) for residents on the 
upper floor was slightly higher than that (16 to 49 
PFU) on the lower floor. Fi et alS 21 ' explained the 
reason as that viral particles were probably trans¬ 
mitted with upflow air movements between flats in 
this block. He and his colleagues* 23 ' have confirmed 
through air modeling that the exhaust fans propelled 
virus-laden air into an outside air shaft, where it was 
carried upward by natural air currents and into other 
flats on the upper floor through open windows. Based 
on our estimate of dose, this transmission mode ap¬ 
pears highly plausible, although we have been unable 
to locate data on viral load in either fecal matter or 
wastewater. 

4. DISCUSSION 

4.1. Validation of Developed Model 

Although we believe that there are no additional 
published data sets available for model development 
for SARS-CoV other than those in Table I, partial 
information supporting our analysis has been pro¬ 
vided by some researchers. Roberts et alS 29 ^ devel¬ 
oped mouse-adapted SARS-CoV through 15 serial 
passages in the respiratory tract of young BAFB/c 
mice and they observed the lethality of mice depen¬ 
dent on the dose of the virus. Mice receiving a dose 
higher than 10 3 9 TCID 50 of the virus died and the 
50% lethal dose (FD 50 ) was 10 46 TCID 50 . In con¬ 
trast, ID 50 of 2.8 x 10 2 PFU estimated with the devel¬ 
oped model is low. Since ordinary SARS-CoV can¬ 
not cause the death of mice, there is a possibility that 
the virulence of the virus was still low even after its 
adaptation to mice. 

Nagata et alS 30) reported that cynomolgus mon¬ 
keys inoculated intranasally with SARS-CoV of 10 6 
TCID 50 did not show any clinical sign and symp¬ 
toms, while slight histopathological changes and virus 
antigen-positive cells were detected. Kuiken et a/.* 31 ' 
and McAuliffe et alS 32 ' 1 also demonstrated the pul¬ 
monary replication of SARS-CoV without any se¬ 
vere symptoms in monkeys receiving intranasal in¬ 
oculation of 10 6 TCID 50 viruses. The fact that as 
many as 10 6 viruses did not cause lethality in mon¬ 
keys means that the susceptibility to this virus may 
be dependent on host species even within primates. 
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Table III. Estimated Dose of SARS-CoV for Residents in Amoy Gardens During the Outbreak in 2003 


(a) Estimation Based on the Assumed Number of Residents 


Floor (Level) 

Number of Reported Cases 

Number of Residents* 

Attack Rate 

Estimated Dose (PFU) 

Lower (4-13) 

12 

320 

0.038 

16 

Middle (14-23) 

46 

320 

0.144 

63 

Upper (24-36) 

41 

416 

0.099 

42 

Overall 

99 

1056 

0.094 

40 

(b) Estimation Based on the Number of Affected Flats 



Floor (Level) 

Number of Affected Flats” 

Total Number of Flats 

Attack Rate 

Estimated Dose (PFU) 

Lower (4-13) 

9 

80 

0.113 

49 

Middle (14-23) 

26 

80 

0.325 

160 

Upper (24-36) 

26 

104 

0.250 

117 

Overall 

61 

264 

0.231 

107 


’Assumed that each flat has four residents. 

’’Number of flats where at least one case was reported. 


Of course, this model originated from data sets 
for mice not for humans. As Mizgerd and Skerred/ 33 ’ 
pointed out, there are the differences between mice 
and humans relevant to pneumonia. However, since 
it is not realistic to obtain the data sets for human in¬ 
fection with SARS-CoV by challenging to volunteers 
due to its high mortality, this model is still a valuable 
method to evaluate the risk of SARS. In prior work, 
it has been shown that data from animal experiments 
provide reasonable estimates for human susceptibil¬ 
ity in inhaled pathogens/ 34,35 ) 

4.2 New Findings About SARS Epidemic Through 

Model Applications 

It is widely considered that SARS is a respira¬ 
tory illness. However, SARS patients may also ex¬ 
hibit gastrointestinal symptoms, splenic atrophy, and 
lymphadenopathy/ 7 ’ Among them, diarrhea is a very 
frequent finding in SARS patients (30% to 40% of 
patients)/ 7 ’ Cheng et alS 36) reported that geometric 
mean titer for SARS-CoV in feces of patients ranged 
from 5.1 x 10 1 copies/mL over 23 days after the on¬ 
set of illness (dpi) to 8.9 x 10 4 copies/mL between 
12 and 14 dpi. Number of gene copies has not been 
related yet to the plaque-forming capability of this 
virus; however, there is no doubt that a huge amount 
of infectious viruses was excreted from patients for 
a certain period. According to the DOH report, at 
the beginning of the outbreak in Block E, a part 
of infectious viruses excreted from the index case 
must have reached other households through the 
waste pipe and ventilation systems. Considering this 
virus transportation, we can imagine that residents 


at the households inhaled small doses of SARS-CoV 
(Table III) as estimated here by model applica¬ 
tion, although there are several unknown factors 
such as virus dispersion and inactivation during 
transportation. 

4.3. Significance of Proposed Model for Preventing 

Reemerging SARS 

The SARS pandemic was brought under con¬ 
trol through a concerted global effort, and by July 
5, 2003, no further human-to-human transmission 
took place. However, there are still several possibil¬ 
ities that might lead to the reemergence of SARS 
in humans/ 37 ) The potential source of SARS-CoV 
might come from infected animals circulating in 
the geographical region, as highlighted by the four 
community-acquired SARS cases between Decem¬ 
ber 2003 and January 2004 in Guangdong, China. In 
these recent cases, fortunately, all patients only de¬ 
veloped mild symptoms and secondary transmission 
did not occur since the animal virus had not fully 
adapted in humans yet. In addition to SARS-CoV, 
two novel Co Vs (HCoV-NL63 and HCoV-HKUl) 
that cause respiratory illness in humans were iden¬ 
tified after the SARS outbreak. These facts indicate 
the possibility that a pandemic of SARS and other 
coronavirus infection may be brought by rapid viral 
evolution. 

5. CONCLUSIONS 

We proposed the exponential model with k — 4. 
1 x 10 2 as a dose-response model for SARS coro¬ 
navirus (SARS-CoV) based on the available data 
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present and would help us predict the reemergence 
of SARS in the future. 
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APPENDIX: DATA SETS FOR MODEL DEVELOPMENT 


No. 

Virus 

Host Animal 

Host Age 

Dose 

Tested 

Positive 

Negative 

1 

rSARS-CoV 

tgMice 

? 

240 

3 

1 

2 





800 

3 

3 

0 





2,400 

2 

2 

0 





12,000 

6 

6 

0 

2 

MHV-1 

Mice 

6 to 8 weeks 

5 

5 

0 

5 





50 

5 

1 

4 





500 

5 

3 

2 





5,000 

5 

5 

0 

3 

HCoV-229E 

Humans 

18 to 50 years 

4 

5 

2 

3 





5 

6 

1 

5 





16 

9 

6 

3 





31 

6 

4 

2 

4 

MHV-S 

Mice 

3 days 

20 

4 

0 

4 





200 

4 

1 

3 





2,000 

4 

4 

0 





20,000 

4 

4 

0 





200,000 

4 

4 

0 

5 

MHV-S 

Mice 

1 week 

200 

4 

0 

4 





2,000 

4 

1 

3 





20,000 

4 

4 

0 





20,0000 

4 

4 

0 

6 

MHV-S 

Mice 

2 weeks 

2,000 

4 

0 

4 





20,000 

4 

0 

4 





200,000 

4 

1 

3 

7 

MHV-2 

Mice 

3 weeks 

700 

5 

0 

5 





7,000 

5 

0 

5 





70,000 

5 

0 

5 

8 

MHV-2 

Mice 

4 weeks 

6,000 

5 

0 

5 





60,000 

4 

0 

4 





600,000 

4 

0 

4 


sets. With this model, the doses of SARS-CoV cor¬ 
responding to 10% and 50% responses (illness) were 
estimated at 43 PFU (95% Cl = 20 to 81 PFU) and 
280 PFU (95% Cl = 130 to 530 PFU), respectively. 
The estimated infectivity of SARS-CoV was compa¬ 
rable to those of HCoV-229E, known as an agent 
of human common cold, and of some animal coro- 
naviruses (MHV-S and HEV-67N) belonging to the 
same genetic group as SARS-CoV. 

The developed model was applied to the anal¬ 
ysis of the epidemiological data of the SARS out¬ 
break that occurred at an apartment complex in 
Hong Kong in 2003. From the reported number of 
cases, it was revealed that the apartment residents 
would be exposed to a dose of SARS-CoV between 
16 and 160 PFU per person, which depends on the 
floor, during the outbreak. 

Although the susceptibility to SARS-CoV seems 
to be host-dependent, the developed model is the 
sole dose-response model for SARS-CoV at the 
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No. 

Virus 

Host Animal 

Host Age 

Dose 

Tested 

Positive 

Negative 

9 

HEV-67N 

Mice 

1 week 

10 

5 

3 

2 





100 

5 

5 

0 





1,000 

5 

5 

0 





10,000 

5 

5 

0 

10 

HEV-67N 

Mice 

4 weeks 

10 

5 

0 

5 





100 

5 

0 

5 





1,000 

5 

2 

3 





10,000 

5 

5 

0 





100,000 

5 

5 

0 

11 

HEV-67N 

Mice 

8 weeks 

10 

5 

0 

5 





100 

5 

0 

5 





1,000 

5 

2 

3 





10,000 

5 

5 

0 





100,000 

5 

5 

0 

12 

HEV-67N 

Rats 

1 week 

20 

5 

3 

2 





200 

5 

5 

0 





2,000 

5 

5 

0 





200,000 

5 

5 

0 

13 

HEV-67N 

Rats 

4 weeks 

2,000 

5 

1 

4 





20,000 

5 

5 

0 





200,000 

5 

5 

0 

14 

HEV-67N 

Rats 

8 weeks 

2,000 

5 

3 

2 





20,000 

5 

3 

2 





200,000 

5 

5 

0 

15 

IBVA-5968 

Chicks 

9 weeks 

0.32 

11 

0 

11 





3.2 

11 

1 

10 





32 

11 

1 

10 





320 

11 

0 

11 





3,200 

11 

1 

10 





32,000 

11 

1 

10 
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